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LONG  TERM  GOALS 

The  long  term  goal  of  this  project  is  to  develop  efficient  inversion  algorithms  for  successful 
geoacoustic  parameter  estimation  and  source  localization,  exploiting  (fully  or  partially)  the  physics  of 
the  propagation  medium.  Algorithms  are  designed  for  geoacoustic  inversion  via  the  extraction  features 
of  the  acoustic  field. 

OBJECTIVES 

•  Achieve  accurate  and  computationally  efficient  geoacoustic  inversion  and  source  localization  by 
designing  estimation  schemes  that  combine  acoustic  field  modeling  and  statistical  modeling. 

•  Develop  methods  for  passive  localization  and  inversion  of  environmental  parameters  that  select 
features  of  propagation  that  are  essential  to  model  for  accurate  inversion. 

•  Implement  sequential  filtering  methods  that  provide  dynamic  and  efficient  solutions  for  the  first 
two  objectives. 

APPROACH 

Continuing  efforts  from  previous  years,  we  worked  with  Bayesian  approaches  applied  to  sound  signals 
for  the  extraction  of  acoustic  features  using  a  combination  of  physics  and  statistical  signal  processing. 

One  of  the  topics  approached  this  past  year  was  geoacoustic  inversion  using  arrival  time  estimation  for 
propagation  in  multipath  environments  with  sequential  Monte  Carlo  methods,  namely,  particle 
filtering.  Some  of  our  results  have  been  presented  in  [1,2,3].  The  idea  is  to  estimate  accurately  arrival 
times  and  corresponding  amplitudes  of  sound  paths  in  shallow  water  environments.  Arrival  times  as 
well  as  amplitudes  of  corresponding  paths  provide  critical  information  on  the  geometry  and 
environment  in  which  sound  propagates.  Even  small  discrepancies  in  arrival  time  estimates  can  lead 
to  significant  errors  to  inferences  made  on  inversion.  Accurate  estimation,  however,  linked  with 
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physical  models  for  sound  propagation  in  the  ocean,  can  lead  to  accurate  geoacoustic  inversion.  In 
addition  to  simple  point  estimates,  our  method  facilitates  the  computation  of  posterior  probability 
density  functions  (PDFs)  for  arrival  times  which  can  be  used  for  uncertainty  estimation  in  geoacoustic 
inversion. 

More  work  was  also  done  on  identifying  arrival  times  and  amplitudes  of  distinct  frequencies  (within  a 
single  mode  or  across  different  modes).  Those  characteristics  provide  significant  information  on 
properties  of  the  propagation  medium  and  source  location.  The  importance  of  modal  arrival  times  and 
amplitudes  in  geoacoustic  inversion  and  source  localization  has  been  discussed  in  [4,  5,  6,  7].  Towards 
this  goal,  the  method  developed  by  the  PI  for  modal  decomposition  of  a  received  signal  and  modal 
arrival  time  estimation  [8,9]  was  applied  for  geoacoustic  inversion. 

WORK  COMPLETED 

The  described  approaches  were  improved  to  provide  more  accurate  estimates  of  arrival  times  and 
geoacoustic  parameter  estimates.  Specifically,  in  addition  to  conventional,  forward  particle  filters, 
smoothing  backward  mechanisms  were  employed,  reducing  uncertainty.  The  arrival  time  tracking 
method  was  applied  to  synthetic  and  real  data  both  for  source  localization  and  sound  speed  estimation 
in  a  seabed  layer.  The  method  was  also  adapted  for  use  with  passive  fathometer  data  [10,1 1],  tracking 
seabed  reflectors.  The  particle  filtering  method  for  tracking  trajectories  in  time-frequency 
representations  was  combined  with  an  acoustic  model  for  sound  speed  inversion. 

RESULTS 

Results  that  we  have  previously  reported  illustrated  how  we  can  extract  arrival  time  estimates  from 
sound  signals  in  the  ocean  that  are  more  accurate  than  those  computed  with  conventional  methods. 
However,  perfonning  geoacoustic  inversion  with  arrival  time  estimates  obtained  with  our  method  for 
some  cases,  we  noticed  that  uncertainty  in  geoacoustic  estimate  results  can  still  be  significant.  An 
example  is  shown  in  Fig.  1,  where  using  multipath  arrivals  extracted  with  particle  filtering,  we 
estimated  sound  speed  in  a  seabed  layer.  Although  the  mode  of  the  PDF  is  very  close  to  the  true  sound 
speed  value  (1645  m/s),  substantial  uncertainty  is  present,  manifested  as  an  extensive  spread  of  the 
PDF. 
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Figure  1:  PDF  of  sediment  sound  speed  estimated  using  arrival  times  of  ray 
paths  extracted  with  particle  filters. 


To  improve  our  geoacoustic  parameter  estimates,  we  further  extended  our  arrival  time  tracking  method 
to  obtain  tighter  PDFs  on  the  arrival  times.  Reducing  the  variance  in  the  arrival  time  PDFs  is  directly 
related  to  limiting  the  variance  in  environmental  parameter  estimates.  We  observed  that  a  portion  of 
the  uncertainty  in  our  arrival  time  estimation  is  a  result  of  limited  prior  information  in  the  initial  states 
of  the  filtering  process.  To  mitigate  this  problem,  we  implemented  backward-moving  filters,  that 
improved  results  on  initial  states  by  moving  backwards  after  estimation  at  all  states/receivers  was 
complete  [12].  Figure  2(a)  shows  the  arrival  time  PDFs  for  a  specific  ray  path  at  a  receiving  phone. 
The  correct  arrival  time  is  49.  The  PDF  that  we  calculate  is  multimodal  with  a  large  variance.  The 
main  mode  is  close  to  the  true  arrival  time.  Figure  2(b)  shows  the  PDF  for  the  same  path  after 
smoothing.  Multimodality  in  the  PDF  has  been  removed  and  variance  has  been  reduced  (however, 
there  seems  to  be  a  small  trade-off  between  uncertainty  reduction  and  a  bias,  which  shifts  the  main 
mode  slightly  to  the  right). 
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Figure  2:  PDFs  of  arrival  time  estimates  of  the  same  path  for  (a)  a  simple  particle  filter 
and  (b)  a  particle  filter  combined  with  smoothing. 

Similar  to  the  example  of  Figure  2  but  now  using  arrival  time  PDFs  obtained  after  smoothing,  we 
calculated  new  estimates  of  the  sound  speed  in  the  layer.  The  new  PDF  is  shown  in  Figure  3.  The 
probability  density  is  now  concentrated  around  the  true  sound  speed  value  with  a  much  reduced 
variance. 


Figure  3:  PDF  of  sediment  sound  speed  estimated  using  arrival  times  of  ray  paths 
extracted  with  particle  filters  after  smoothing. 
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Further  work  was  carried  out  by  applying  our  method  to  real  data.  Specifically,  we  worked  with  data 
from  the  SW06  experiment,  collected  at  the  16-element  MPL-VLA1  array  (because  of  low  SNR  at  the 
15th  and  16th  phones,  data  at  only  the  14  lower  phones  were  used)  [2],  The  source  signal  was  a  linear 
frequency  modulated  pulse  with  frequencies  between  100  and  900  Hz.  Figure  4  (a)  illustrates  ray  path 
arrival  time  PDFs  in  a  case  where  the  number  of  paths  arriving  within  a  certain  observation  time 
window  is  varying.  Figures  4  (b),  (c),  (d)  show  PDFs  for  the  number  of  observed  arrivals  within  the 
given  window  at  three  different  phones. 


Number  of  arrivals 


Time  (ms) 


Number  of  arrivals 


Figure  4:  (a)  PDFs  for  arrival  times  of  three  paths  calculated  with  particle  filtering  applied 
to  the  time  series  of  Fig.  1.  (b)  Probability  Mass  Functions  for  the  number  of  arrivals 
in  an  observation  window  for  receivers  9,  8,  and  7  [2], 


As  previously  mentioned,  this  approach  is  highly  applicable  to  data  produced  by  passive  fathometer 
data  processing  [10,1 1].  The  presence  of  wavelets  in  the  estimate  of  the  medium’s  Green’s  function 
corresponds  to  the  location  of  reflectors  in  the  seabed.  Figure  5  shows  how  the  particle  filtering 
method  tracks  reflectors  from  data  from  the  Boundary  2003  Experiment,  after  both  conventional  and 
MVDR  processing.  This  work  was  done  in  collaboration  with  Drs.  Peter  Gerstoft  and  Caglar  Yardim 
from  Scripps. 
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Figure  5:  (a)  Conventional  and  (b)  MVDR  fathometer  results  for  Boundary  2003  data,  (c) 
Conventional  and  (d)  MVDR  fathometer  results  for  Boundary  2003  data  with  particle-filtering 

tracks  superimposed  (white  dots). 

Inversion  was  also  conducted  using  time-frequency  particle  filtering  [8].  Figures  6  and  7  show  sound 
speed  PDF  estimates  (at  one  depth  in  the  water  column  and  within  the  first  sediment)  for  synthetic 
data.  The  ellipses  indicate  the  true  values  of  the  sound  speed.  The  inversion  was  done  using 
frequency-arrival  time  estimates  and  linearization.  Ill-conditioning  issues  presented  us  with  challenges 
which  we  are  currently  investigating. 


Figure  6:  PDF  of  water  column  sound  speed  at  a  selected  depth  calculated  using 
particle  filtering  and  time-frequency  analysis. 
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Figure  7:  PDF  of  sediment  sound  speed  calculated  using  particle 
filtering  and  time-frequency  analysis. 


IMPACT 

The  significance  of  accurate  arrival  time  estimation  in  geoacoustic  inversion  has  been  studied  with 
several  methods  designed  for  producing  geoacoustic  parameter  estimates  and  measures  of  the 
uncertainty  in  the  estimation  process.  The  reliability  of  these  methods  is  intimately  tied  to  the  ability  to 
accurately  extracting  and  identifying  arrival  times.  The  new  methods  facilitate  this  extraction  and  the 
association  between  paths/modes  and  detected  arrivals  and  also  produce  posterior  probability 
distributions  of  modal  frequencies  or  arrival  times.  These  distributions  can  then  be  employed  to 
quantify  uncertainty  in  the  estimation  of  geoacoustic  parameters  as  demonstrated  in  Figures  3,  6,  and  7. 
The  success  of  the  developed  approach  in  reflector  tracking  using  passive  fathometer  data  is 
encouraging,  as  this  is  an  emerging  field  with  a  vast  potential  for  the  understanding  of  the  seabed 
structure. 

RELATED  PROJECTS 

The  PI  is  collaborating  with  Drs.  Yardim  and  Gerstoft  on  sequential  filtering  in  ocean  acoustics.  A 
collaboration  is  also  underway  with  Dr.  Lisa  Zurk  and  Dr.  Helen  Ou  (Portland  State  University)  on 
tracking  as  it  relates  to  the  waveguide  invariance  principle.  Additionally,  the  PI  is  collaborating  with 
Dr.  Leon  Cohen  on  comparing  numerical  and  analytical  descriptions  of  dispersion  in  ocean 
environments. 
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